Impurity Quantum Phase Transition in a Current-Carrying <i-Wave Superconductor 
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We study an Anderson impurity embedded in a d-wave superconductor carrying a non-zero mo- 
mentum of Cooper pairs driven by magnetic flux through the Aharonov-Bohm effect. The zero- 
temperature impurity behavior is investigated by using the numerical renormalization group method 
in calculating the energy flows, the impurity spectral functions, and the effective magnetic moments. 
These results explicitly show that the local impurity state can be changed dramatically from the 
asymptotically free spin to the completely screened one by tuning the magnetic flux. 

PACS numbers: 71.10.Hf, 71.27. +a, 75.20.Hr, 71.28.+d 



The behavior of a single magnetic impurity in corre- 
lated electron systems has attracted intensive interest 
in condensed matter physics While the problem is 
well studied when the host is a simple metal 0, the 
correlations among the host electrons may lead to much 
complicated impurity behavior deviating from the Fermi 
liquid properties Specifically, the Kondo screened 
state, i.e., an entanglement state between the impurity 
spin and a conduction electron, may become unstable 
at zero temperature upon the depletion of the density of 
states (DOS) at the Fermi levels, or the change of certain 
nonthermal parameters. Such kind of impurity quantum 
phase transition (IQPT) may take place in metallic sys- 
tems with a pscudo gap The response of a known 
impurity behavior upon the change of effective couplings 
could provide important information of the host itself, 
and thus can be used to probe the ground state and 
low-energy physics of the host electrons On the 

other hand, within a known host, the different response 
of a quantum impurity with internal dynamical degrees 
of freedom versus a static impurity, to an external con- 
trol parameter, may shed insight on the role played by 
a doped impurity in a correlated electron medium [B|, [y] , 
which is a relatively rarely explored area. In this Let- 
ter, we study the property of an Anderson impurity em- 
bedded in a d-wavc superconductor. Our analysis based 
on the numerical renormalization group (NRG) method 
suggests that an IQPT can be induced in a d-wave super- 
conductor carrying persist currents in a controllable way. 
The obtained results can be tested by scanning tunnel- 
ing microscopy (STM), which has been used to measure 
the local electron structure around a doped impurity in 
conventional Q and unconventional superconductors Q 
in the absence of a current flow. 

The d-wave superconductor, usually being of quasi 
two-dimensional layered structure, could be designed 
with a geometry of hollow cylinder, as shown schemat- 
ically in Fig. [TJ The persistent (super-)current can be 
generated by piecing a magnetic flux through the axial 
of the cylinder, due to the Aharonov-Bohm effect 0. 








FIG. 1: (Color) An Anderson impurity embedded in a two- 
dimensional d-wave superconductor made in the hollow cylin- 
der geometry. A persistent current q is generated by the 
Aharonov-Bohm effect which can be tuned by varying the 
magnetic flux $ passing through the cylinder. The local elec- 
tron structure around the impurity can be measured by scan- 
ning tunneling microscopy. 

We can model the above impurity problem by coupling 
an Anderson impurity to the d-wave superconductor car- 
rying a current, using the following Hamiltonian 

~H = WbCS + Wimp + Whybridj (1) 

where Hbcs = E k ,a 6c+q4 CT c k<r + EkI A kCk T cL ki + 
h.c], Wi m p = Y, a e dd\d a + UndfUdi, and ^hybrid = 
"yjvT Ek crl^kC^rfo- + h.c.]. Here Ck<j annihilates one con- 
duction electron of momentum k and spin projection 
cr, while d a annihilates one localized d-electron of spin 
projection a. In the tight-binding approximation, the 
conduction electrons have the normal and d-wave super- 
conducting gap dispersions, £k = — 2t(cosfc x -l-cosfcy) — 
At' cos k x cos k y — ^ and Ak = (Ao/2)(cosfca; — cosk y ), 
respectively. Notice that the energy is measured with 
respect to the Fermi energy. The center-of-mass momen- 
tum, 2q, of a Cooper pair determines the current flow. 
Other parameters e^, U, and Vk are the localized level, 
the on-site Coulomb interaction, and the impurity cou- 
pling, respectively. Nr, is the number of lattice sites. 

Distinct from the case of either a wide-band normal 
metal with a constant density of states, where the im- 
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purity spin being in a strong coupling (SC) limit, or 
a conventional s-wave superconductor with a hard-gap 
everywhere on the Fermi surface, where the impurity 
spin being in the local moment (LM) limit, the exis- 
tence of nodal zero-energy quasiparticles in a d-wave su- 
perconductor 0, have a non-trivial implication to 
the fate of quantum impurity states. Earlier studies by 
taking the DOS with a soft-gap, p(u) ~ \u\ r (r = 1 
for the d-wave superconductor), have shown the exis- 
tence of a critical coupling, separating the LM and SC 
phases [1, 11, H, l]| 0, EH- More intriguingly, we no- 
tice that the low-energy excitations in a d-wave super- 
conductor can readily be proliferated by pumping in a 
supcrcurrent, which should have a significant control of 
the impurity state. The purpose of the present work is to 
demonstrate that the IQPT can be realized by gradually 
increasing the supercurrent. 

In order to solve the problem, we generalize the NRG 
method 0, [l(| to study the impurity properties with 
an arbitrary form of the DOS, which in the present 
case is found by diagonalizing 7Ybcs to be p(E) = 
E k [K q | 2 <5(£ - 4?q) + k, q | 2 <5(£ - E^)], where 
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with Z kiq = (£ k+q - 
q )/2. The electron- and 
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£k- q )/2 and Q kjq = (£ k +q+£k 
hole components of the Bogoliubov-de Gennes eigenfunc- 
tion are given by |w(w)k, q | 2 = [1 ± Qk,q/-E k ° q ]/2. 

In the numerical calculations, we take t=l,t'= —0.2, 
p = —0.78, Ao = 1. The energy is measured in units of 
t = 1 unless specified otherwise. Without loss of gener- 
ality, we take q y = 0. This choice of q = {q x ,q y ) corre- 
sponds to the circular current perpendicular to the axis 
as shown in Fig. [TJ Its magnitude is then related to the 
magnetic flux via q x = 27ry> / '(Nl$o) with <f>o ~ hc/e. In 
Fig. [2l we show the DOS as a function of energy for var- 
ious values of q x . A small intrinsic lifetime broadening 
2 x 10 -4 , and the lattice sites of 2 14 x 2 14 are chosen. In 
the absence of the current, the calculated DOS vanishes 
linearly as expected. However, it increases around the 
Fermi energy as a response to the non-zero current. 

It is sufficient for our purpose to take into account 
only the coupling between the particle-excitations and 
the impurity spin. Indeed, it is argued that the local 
quasi-particlc DOS of the superconductor is the only nec- 
essary ingredient in a number of cases, in particular the 
unconventional superconductors [It], EH ■ For simplicity, 
we therefore ignore the anomalous part and just study a 
modified Anderson model with the impurity coupled to 
the electron-like excitation spectrum. As such, we use 
the following Hamiltonian to the derivation of the NRG 
equations 
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FIG. 2: (Color) Density of states as a function of energy for 
various values of drift momentum q x . 



where we introduced a one-dimensional energy represen- 
tation for the particle-like excitations a\ a with the scaled 
energy e and the band-cut-offs at ±D. The quantities 
g(e) and h(e) are the energy dispersion and hybridiza- 
tion, seli-consistently defined by p(E) and respectively 
as in [3]. 

To handle the arbitrary p(ui), we follow the generic 
scheme of Ref. [lj| to discretize the Hamiltonian given 
by Eq. ((2]), and then solve it by the NRG for different 
values of U and A k at the symmetric point €d = —{7/2, 
with the RG parameter A — 2.5. At each iteration step 
we keep sa 200 ~ 1200 states, depending on the quantities 
calculated, and the total number of iterations is iV max w 
100. Throughout the work, a value of D = 2 is chosen. In 
the following, we briefly discuss the results of the energy 
flows, the impurity density spectra, as well as the effective 
impurity moment in various cases. 

Actually even in the absence of the current, there 
should be an IQPT from the SC state to the LM state 
upon increasing repulsive U at the finite impurity cou- 
pling. Fig. [3] shows the low lying energy levels as a func- 
tion of odd number of iteration N for various values of 
U in increasing order and for fixed (and k-independent) 
V = 0.01. For small values of U (see, e.g., Panel (a)), 
for small TV, the levels start and remain close to the free 
orbital fixed point. For larger N, they make a rapid 
crossover to the LM fixed point. Due to the marginal 
nature of the (i-wave superconduting host, they remain 
near the LM fixed point for a finite number of iterations 
as manifested by a plateau, and then gradually crossover 
to the SC limit. This length of the plateau increases with 
increasing U (see Panels (a),(b)). For the large values U, 
the low lying energy levels always flow to a LM fixed point 
(see Panels (c),(d)), suggesting the impurity state is in 
the localized regime. Correspondingly, we show in Fig. [4] 
the impurity spectral function and the effective magnetic 
moment for fixed V = 0.01. As can be seen in Fig.dJ the 
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FIG. 3: (Color) Flow diagrams for the low-energy many- 
body excitations obtained from the numerical renormaliza- 
tion group for various Coulomb U with drift momentum 
q x = q y = 0, N (odd) is the number of iterations of the 
NRG procedure, A is the discretization parameter. Solid lines: 
(Q,S) = (1,0), dashed lines: (Q,S) = (2,1/2), with Q and 
S being the quantum numbers of the charge and spin respec- 
tively. 
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FIG. 4: (Color) The impurity spectral function at T — 
(main panel) and the temperature dependence of the effective 
magnetic moment (inset) for various Coulomb U, s/ = —(7/2 
with drift momentum q x = q y =0. 



spectral density exhibits a narrow Kondo resonance peak 
at the Fermi energy in the SC state while shows only a 
gap-like behavior in the LM state. Also from the inset 
of Fig. [U one can see that the screened moment, which 
is essentially zero in the SC state, saturates at the maxi- 
mal value 0.25 as U approaches to a larger value. These 
results provide a base to study the effect of non-zero q 
starting from the LM phase. 

We then calculate the similar quantities for fixed U = 
0.1, and V = 0.01 with q y = and various q x . As shown 
in Fig.0 the length of plateau in the energy flow becomes 
shortened with the increasing q x (i.e., the supercurrcnt 



FIG. 5: (Color) Flow diagrams for the low-energy many- 
body excitations obtained from the numerical renormalization 
group for Coulomb U = 0.1 with various drift momentum q x , 
N (odd) is the number of iterations of the NRG procedure, 
A is the discretization parameter. Solid lines: (Q, S) = (1, 0), 
dashed lines: (Q,S) = (2, 1/2). 



flow). This evolution of energy flows from the LM line- 
shape to the SC one is similar to the soft-gap Anderson 
impurity model. It indicates that the impurity state orig- 
inally staying in the LM regime can be driven into the 
SC state. In the present case, the transition takes place 
around q x w 0.003, though it is not so sharp as com- 
pared with the soft-gap model. The impurity spectral 
function begins to show a central peak at zero energy as 
we increase q x across this critical value (see Fig. |6]) . This 
feature is similar to the soft-gap Anderson model where 
the impurity spectral function diverges at zero energy 
for the SC and quantum critical phases (Td| . Remark- 
ably, the sum rule is within 95% accuracy in our case, 
which is better than that in the soft-gap model. The 
IQPT is also indicated in the effective moment, plotted 
in the inset of Fig. [6] It shows that the effective moment 
approaches to the free moment value 0.25 for q x < 0.003, 
and approaches to zero for q x > 0.003. The later case in- 
dicates a SC phase where the impurity spin is completely 
screened. Our numerics also shows that the increase in 
the number of states, kept in the iterations, improves the 
proxy of the saturated values of the effective moment to 
0.25 (0) in the LM (SC) limit. 

Our NRG results demonstrate the existence of finite 
critical point around q x>c « 0.003 when q y = 0, U = 0.1, 
and V — 0.01. The impurity is in the LM and SC phases 
for q x < q x<c and q x > q x ^ c , respectively. We find that 
this conclusion is robust if, for q x = q y = 0, V varies in 
the LM phase. In general, q X)C will increase with U for 
fixed V. The precise value of q x ^ c requires extremely high 
precision of the DOS of the low-energy excitations, which 
could be achieved within 10 -16 ~ 10~ 12 for the soft-gap 
Anderson model. However, for the present model, such 
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will emerge very close to the Fermi energy while those in 
the large-gap domains are still in the LM with no Kondo 
resonance. These interesting phenomena, which can be 
observed by STM, should serve as a direct test of our 
prediction. 
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FIG. 6: (Color) The impurity spectral function at T — 
(main panel) and the temperature dependence of the effective 
magnetic moment (inset) for Coulomb U = 0.1, £/ = —{7/2 
with various drift momentum q x . 



high accuracy of the DOS needs sufficiently small intrin- 
sic lifetime broadening and sufficiently large system size 
Nl, so that a detailed analysis on the scaling behavior 
in the critical phase is challenging. Nevertheless, it is 
sufficient for our purpose to predict a critical point sep- 
arating the SC and LM phases by tuning the current or 
the magnetic flux. 

Our findings have several implications: (i) In a d-wave 
superconductor, quasiparticle resonances can be induced 
around a doped static impurity in the strong scatter- 
ing limit. Upon the flowing of a supercurrent, the res- 
onance peak is suppressed in amplitude and broadened 
in width [2fj| . In contrast, for the quantum impurity, 
a Kondo resonance emerges when it is driven into the 
SC phase by the supercurrent. This distinction of the 
response to the supercurrent can help us to decipher 
whether a doped atom plays the role of a static or quan- 
tum impurity in a d-wave superconductor, (ii) In high-T c 
cupratcs like Bi2Sr2CaCu20s+5, the nanoscale inhomo- 
geneity with the existence of small- and large-gap do- 



2 ll . l22j has been ubiquitously observed by the 



STM. When a magnetic impurity is doped into the sys- 
tem, it is in the LM phase in the absence of a super- 
current. In this case, the impurity plays the role of a 
weak potential scatter [23|, and the resonance as due to 
the quasiparticle scattering is located far away from the 
Fermi energy. This resonance should be observed in both 
types of domains. When a supercurrent flows in the sam- 
ple, the impurities in the small-gap domains should be 
first driven into the SC phase and a Kondo resonance 
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